Equivalence and limitations of error suppression techniques 
for adiabatic quantum computing 
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While adiabatic quantum computation (AQC) possesses some intrinsic robustness to noise, it is ex- 
pected that a form of error control will be necessary for large scale computations. Error control ideas 
developed for circuit-model quantum computation do not transfer easily to the AQC model and to 
date there have been two main proposals to suppress errors during an AQC implementation: energy 
gap protection and dynamical decoupling. Here we show that these two methods are essentially 
equivalent and can be analyzed within the same formalism. We analyze the effectiveness of such 
error suppression techniques and identify critical constraints on the performance of error suppression 
in AQC, suggesting that error suppression by itself is insufficient for fault-tolerant, large-scale AQC 
and that a form of error correction is needed. 



Adiabatic quantum computation (AQC) [1, 2] is often 
praised for its inherent robustness to both dephasing 
and energy relaxation, phenomena known to plague the 
majority of quantum computing paradigms. Numerous 
studies, however, have shown that single qubit noise is 
capable of driving undesirable transitions from the adi- 
abatic ground state; e.g. Refs. [3-7]. This observa- 
tion prompted the formulation of three error suppression 
techniques, each leveraging the error detecting properties 
of quantum stabilizer codes [8]: energy gap protection 
(EGP)[9], in which the addition of the stabilizer gen- 
erators to the system Hamiltonian causes errors to incur 
large energetic penalties; dynamical decoupling (DD)[10], 
whereby stabilizer generators are applied periodically as 
unitary operators, rcfocusing errors much like traditional 
spin echos; and Zeno effect suppression [11], which in- 
hibits errors by frequent measurement of the stabilizer 
generators. These three techniques apparently operate 
by very different physical mechanisms. However, Fac- 
chi et al. [12] have shown that both Zeno suppression 
and DD may be viewed as limiting cases of a more gen- 
eral mathematical framework. In this work, we show the 
DD and EGP may themselves be unified under a single 
formalism, and that these two methods are effectively 
equivalent in their error suppression power. 

In conventional quantum computing models (e.g. the 
circuit-model) it is well understood that such error sup- 
pression techniques by themselves are insufficient for 
fault-tolerant quantum computing. From a thermody- 
namic perspective this is because error suppression alone 
does not provide a mechanism to remove the entropy gen- 
erated by errors from the encoded system. Since the 
thermodynamic argument is independent of the compu- 
tational model, it is expected that error suppression alone 
is insufficient for fault-tolerant quantum computing in 
the AQC model as well. We provide strong evidence to 
support this expectation using several arguments. 



An outline of the paper is as follows: first we prove 
the equivalence of the EGP and DD error suppression 
techniques for AQC using a formal correspondence be- 
tween the Hamiltonians describing system evolution in 
an interaction picture. Then we show using a simple ar- 
gument that, unlike circuit model quantum computing, 
encoded AQC does not tolerate any excursions outside 
the codespace, and errors in general are unrecoverable. 
This places stringent requirements on error suppression 
and correction timcscales in AQC. Finally, we derive a 
master equation of non-Lindblad form that describes en- 
coded adiabatic evolution by exploiting properties of the 
stabilizer encoding. This master equation makes appar- 
ent the dynamical mechanism behind error suppression 
and also allows us to analyze the scaling limitations of 
error suppression. 

Encoded AQC - Consider a Hamiltonian acting on 



a tensor product Hilbert space, H 



He 



describ- 
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ing a system undergoing adiabatic quantum evolution 
while coupled to an external environment/bath: H(t) = 
#aqc(£) + YJj=i E j ® Bj + H B - Here H A qc acts on the 
system and performs the adiabatic evolution, Bj and Ej 
are bath operators and single qubit Pauli error operators, 
respectively, that define the system-environment inter- 
action. Hb is the bath Hamiltonian. Time-dependent 
eigenstates of Haqc may be labeled as \n,k) according 
to their principle quantum number n and index k distin- 
guishing any degeneracy. 

We assume that the evolution is sufficiently slow that 
we may safely ignore diabatic transitions. All errors are 
then due to the system-bath interaction which is able to 
induce transitions from the adiabatically evolving ground 
state and therefore cause the computation to fail. To sup- 
press these transitions using any of the techniques men- 
tioned above requires encoding the system in an error de- 
tecting stabilizer code [8]. Typically, the code is chosen 
to be the smallest code for which each Ej in the system- 
bath interaction anti-commutes with at least one of the 
stabilizer generators. This encoding will enlarge the sys- 
tem's Hilbert space by a factor of 2 Ng , where N g is the 
number of stabilizer generators of the code. The physical 



2 



operators, a x ,a y , a z in Haqc are replaced by the code's 
logical operators, X, Y, Z, and a time-dependent system 
control Hamiltonian composed of elements from the sta- 
bilizer group (specified later) is added to implement any 
desired error suppression. The encoded Hamiltonian is 
then: 

H{t) = H AQC (t) + H c (t) + Y,E J ®B J +H B ( 1 ) 



We have assumed that the system-bath interaction re- 
mains qualitatively the same after the encoding, but 
is extended to N e > n e terms to correspond the 
larger system size. States of the encoded system 
may now be labeled by the same two quantum num- 
bers as before, but with N g additional quantum num- 
bers given by the eigenvalues of the stabilizer genera- 
tors, S m \n, k; 81,82, . . . , s Ng ) = s rn \n, k; si, s 2 , . . . , s Ng ), 
where S m is a generator of the stabilizer group and 
s m = ±1. We shall refer to the codespace as the set 
of states for which all the stabilizer eigenvalues arc +1. 

For the following discussion it will be useful to define 
U n {t u t 2 ) = exp + H ji* 2 dsH n {s)), for n e {AQC,B,C} 
(+ denotes positive time ordering but this is unnecessary 
except for n = AQC), and U n {t) = U n (t,0). Note that 
[Haqc(s),H c (s')} = Vs, s', because H A qc only con- 
tains logical operators of the code and Hq only contains 
stabilizer terms, and hence these unitaries all commute 
with each other. The following notation is used for an 
operator A in an interaction picture with respect to the 
control: A{t) = ul(t)Ul(t)A U c (t)U B {t), which is typ- 
ically called the toggling frame. Evolution of states in 
this frame is generated by the toggling frame Hamilto- 
nian: H{t) = U ] c {t)U B {t) {H{t) -H c - H B ) U B {t)U c {t). 
Dynamical decoupling - Dynamical decoupling con- 
trols are chosen to sequentially apply the generators of 
the stabilizer group as unitary operators. The stabilizers 
are applied in a particular order, given by the vector n, 
at times given by K it) 6 Z, so that at time t the most 
recent operator applied to the system was S nK(ty The 
unitary operator defining the toggling frame may then 
be written as, W^? D (<) = flj^ ^Vr ^ s a P r °duct of sta- 
bilizer generators, the operator, W^ D (t), is an element 
of the full stabilizer group and therefore commutes with 
.Saqc (*) ■ I n t ne toggling frame the encoded Hamiltonian 
takes the form: 

H OD (t) = ii aqc (t) + J2 E T^) ® 

And because Ej must either commute or anti-commute 
with each member of the stabilizer group [8], we write: 

E^it) = U° D \t)E 3 US D (t) = i-lf (t) E (2) 

where p{t) = if [E 3 ,Ug D (t)] = and p{t) = 1 if 
{E 3 ,U^it)} = 0. An effective DD cycle is one that 



causes pit) to rapidly alternate between +1 and —1, re- 
sulting in the system-environment coupling being modu- 
lated by a rapidly oscillating function of t. The state 
of the combined system-plus-environment at time t is 



given by: 



*(*) 



exp + (-i J Q dsH BB (s)) 



with 



and we set h = 1 through- 
out. Integrating over the modulation factor (— in 
the exponential has the effect of suppressing the system- 
environment coupling over timescales longer than the DD 
inter-pulse period. This intuition can be made more ex- 
plicit by deriving the average Hamiltonian using a Mag- 
nus expansion[ 1 3] whereupon the oscillating terms may 
be seen to cause a reduction in the average system-bath 
coupling [14]. 

Energy gap protection - The EGP approach proceeds 
by setting the control Hamiltonian equal to a sum of sta- 
bilizer generators, Hq GF '{t) = — ctJ2m=i ^m, with a > 0. 
States in the codespace are then eigenstates of He with 
eigenvalue —aN g , but any state outside the codespace is 
subjected to an energy penalty of at least a which is ex- 
pected to reduce transitions into these states. Since He 
is a function of only the stabilizer generators, UQ GP it) 
again commutes with the code's logical operators which 
comprise i?AQc(i)- so we can write the Hamiltonian in 
the toggling frame as: 

H EGP {t) = H AQC {t) + Ef GP it) ® Bj{t) 

Error operators in the EGP toggling frame can be shown 
to take the form: 



Ef GP it) = E j e^ iat ^ s ^= a 



s„ 



-2iatJ2 



{S m ,E, }=0 ' 



'E 



j • 



(3) 



where the sums are taken over all stabilizer generators S m 
that anti-commute with the error operator Ej . To obtain 
this expression we have exploited the following: (i) the 
stabilizer generators commute with each other, and (ii) 
each generator either commutes or anti-commutes with 
the noise operators: S m Ej = ±EjS m . Let Wj be the 
number of generators that anticommutc with Ej. Then 
the action of this toggling frame Hamiltonian on any 

state, I*) = J2 t Wc) ® W) e 'Wsys ® H cm with |$) in 
the codespace is, 



N, 



H 



EGP 



(*)!*) 



^2iuu at 



Bj{t) 1*) 



Thus the coupling term Ej <S$Bj{t) is modulated by a fac- 
tor of e jat . The evolution of states in the codespace is 
described by exponentiating the integral of this Hamilto- 
nian and integrating over the modulation factors above 
has the effect of suppressing the coupling of the system to 
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the environment. Similar to DD, the error term is modu- 
lated by an oscillating function when acting on states in 
the codespace in the interaction picture. In the case of 
EGP the oscillations are smooth and sinusoidal whereas 
the oscillations for (impulsive) DD are square waves in 
time. One may then mimic a decoupling sequence by 
choosing a so that the EGP oscillations match the fre- 
quency of a DD sequence. Numerical studies provide ev- 
idence that in such cases EGP and DD both yield nearly 
identical evolutions of states in the codespace. And since 
error suppression performance is dictated solely by evo- 
lution of states in the codespace (the subspace of cor- 
rectable states is irrelevant for AQC, as shown in the 
next section) we conclude that EGP and DD are effec- 
tively equivalent strategies for error suppression. 

This in fact motivates generalizations of EGP where 
the weight terms, a, are not constant in time or equal 
across the stabilizer generators. Many decoupling se- 
quences vary the time interval between the pulses; 
UDD [15], for example, chooses the pulse arrival times 
cos(nir/2(N + 1)), where N is the total num- 
ber of pulses in time interval [0,T]. To mimic this UDD 
sequence, where the modulation frequency is not con- 
stant in time, we choose a time dependent weight term, 
a(t) = NT/ \/t(T — t). This approach was used, in a 
slightly different context, by the authors of Ref. [16] 
to produce an effective UDD sequence using continu- 
ous controls. More generally, allowing a to vary in 
time allows the strongest identification between the DD 
and EGP approaches, and a unified treatment of both 
as quantum control protocols. For instance, choosing 
ctj(t) = J2i — *i ) /2 applies Sj at time t\ as a uni- 
tary operator (impulsive DD), but in the EPG formalism. 
Furthermore, this approach naturally lends itself to the 
application of optimal control techniques to choose Uj (t) 
to optimally mitigate the system-bath interaction. 

Error suppression in AQC - We pause to empha- 
size a key difference between encoded circuit-model quan- 
tum computation and encoded AQC. Typical analyses 
of the former consider any state in the correctable sub- 
space as uncorrupted since these states can be decoded 
perfectly at the end of the computation. However, this 
is not true for AQC. This can be seen in a number of 
ways, but perhaps most transparently through an exam- 
ple evolution under a stochastic noise model [17]. Con- 
sider unperturbed evolution of the system until time r, 
at which point there is a correctable error Ej. After 
this error the system evolves unperturbed again until the 
end of the AQC. Thus the overall evolution is: \ip(T)} = 
U A qc(T,t)E j U A qc(t,0) |V' c (0)). Note that we are only 
considering evolution in H sys since this is sufficient and 
we are ignoring evolution by the control Hamiltonian 
since it is inconsequential for the following argument. 
This is equivalent to: \ip(T)) = EjU A Qc{T, 0) \ip c (0)}, 
where Ej = U A qc(T, T)EjU AQC (T, r). That is, ideal 

evolution proceeded by an application of the error Ej. 
Since Uaqc (T, t) is a function of the logical operators of 



the code, Ej is in general a many-body operator (con- 
tains a non-identity Pauli on more than two qubits of 
the system). And hence, although Ej is a detectable 
single-body Pauli error, Ej can induce logical errors that 
are undetectable by the code. Therefore we cannot re- 
cover the state at time T with an application of a detec- 
tion and correction operation even though the original 
error at r was detectable. Error correction (or decod- 
ing) at the end of the computation is unlikely to be ef- 
fective; correctable errors will be transformed to logical 
errors by the Hamiltonian over the course of the compu- 
tation. This is analogous to an error during the imple- 
mentation of a non-transversal gate in the circuit model. 
This scrambling of errors by Haqc (t) places constraints 
on pure error suppression techniques in AQC because it 
means that the suppression has to be strong enough to 
keep the state in the code space (and not just in the cor- 
rectable space). This motivates active error during the 
course of AQC evolution, though it also imposes strict 
conditions on the rate at which error correction must be 
performed. Specifically, one must correct on timescales 
fast compared to the rate at which i?AQc(^) moves an 
erred state out of the correctable space. We can estimate 
the correction timescale, r c , by the demanding that it is 
fast compared to the norm of the AQC Hamiltonian, i.e., 

\\H A Qc(t)\\T c < 1. 

Master equation for encoded AQC - We have shown 
that error suppression in AQC using DD or EGP are for- 
mally equivalent as far as codespace population dynamics 
is concerned. Now we formalize this further by deriving 
a master equation describing effective encoded adiabatic 
evolution when the qubits interact with an environment. 

We begin with the full Hamiltonian in Eq. (1), and as- 
sume that the system-environment coupling is weak com- 
pared to the other terms in this Hamiltonian, and that 
the encoding has been chosen such that each Ej is a de- 
tectable error. In addition to the weak coupling approx- 
imation we also employ the first Markov approximation 
[18]. Crucially, we do not employ the second Markov 
approximation [18] in the derivation and therefore the 
reduced dynamics is able to capture the modification of 
system-environment coupling, and hence decoherence, by 
controls such as dynamical decoupling or energy penalty 
terms. This allows us to derive the following master equa- 
tion in the toggling frame [19]: 



dp(t) 
dt 



i[H AQ c(t),p(t)} 



3 = 1 



Cj^Ej^Ejit^pit) 



C*(T)Ej(t)p(t)Zj(t,T) +h.c 



where Ej(t,T) = U A Qc(t,t - t) Ej(t-r) U AQC (t,t - r), 
and C,j{t) = Tr eliV {W B (T,0)BjU B (T,0)B j a ecl } is the 
bath correlation function. We have assumed that the 
environments of all the qubits are uncorrelated and sta- 
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tionary. Correlated environments can be captured by the 
same formalism but for simplicity we will not do so here. 

Now consider the change in the codespace population 
as a result of these dynamics. Let P = f]m = i |(I + S m ) 
be the projector onto the codespace, and Q = I P. 
Then the change in the codespace population is ^f- = 
Tr{P|fP}. We can show that the dynamical equation 
for this code space population is well approximated by 
[191: 



N e 



dr Re < Cj (r)m 3 - (t, r) x 



(4) 



Tr £7 J -PH J -(i,T)Qp(i)Q - EjEjit, r)Pp(i)P 



dP c (t) 
dt 



where Ej(t,T) = e - iTHA ^ c( -^ Ej e iTHA » c( - t \ and m^r) 
is a modulation function that results from the control, 
and captures its influence on dissipation and decoher- 
ence. For the two cases of EGP and DD, this modulation 
function takes the form: 



EGP/ 



,DD 



(t,r) = (-1)p(*)-p(*- t ) 



(5) 



These modulation functions arc analogous to the fil- 
ter functions derived for describing dynamical de- 
coupling for pure dephasing dynamics [20]. In 
the absence of the adiabatic evolution (i.e. when 
-Haqc = 0) the traces in this equation sim- 
plify to Tr^Sj^rJPp^P} = Tr{Pp(t)P}, and 
Tr^PS^rjQp^Q} = Tr {QipO)Qi}, where Q x is 
a projector onto the subspace of Q that contains states 
one error away from the codespace. This simplification 
allows the derivation of a classical master/rate equation 
for the codespace population: 



dP c (t) 



dt 



«AQC=0 



where r^(t) = 2Re{ J Q d,TCj(T)mj(t 7 T)} and P e \(t) is 
the population in the one-error subspace at time t. 

The effect of the control Hamiltonian can be clearly 
seen in this master equation. The rates T^{t) quantify 
the leakage out of the codespace per unit time. In the 
absence of a control Hamiltonian these rates are sim- 
ply proportional to Re{/ dTCj(r)}, a property of the 
environmental fluctuations alone. However the control 
Hamiltonian, in the case of DD and EGP, has the effect 
of modulating the integrand with an oscillating function 
and hence decreasing its amplitude. To illustrate this 
consider a classical approximation of the environment 
(e.g. the Kubo-Andcrson stochastic model) for which 
the correlation function is purely real, and fix it to be 
exponentially decaying [21], C(t) oc e~ 7 * where 7 is the 
inverse correlation time. Consider EGP, where the mod- 
ulation function is sinusoidal, for which, 



rf(t) cx 



2 (7 — 7e 7 * cos[2aWjt] + 2awje 7 * sm[2awjt]) 



In this expression, the consequences of adding the en- 
ergy penalty terms are summarized by the presence of 
the factor 2awj , which increases with the energy penalty 
(a) and the number of that anti-commute with the er- 
ror (wj). The term has two effects: (i) its presence in 
the denominator decreases the overall rate of population 
leakage, (ii) it increases the oscillation frequency of the 
sinusoidal functions in the numerator, thus decreasing 
the magnitude of integrals of r^{t). Therefore, this cal- 
culation explicitly shows how the control Hamiltonian 
decreases population leakage from the codespace. Note 
that it is possible to achieve a suppression of population 
leakage from the codespace that is exponential in the 
number of qubits by adding the whole stabilizer group 
as the penalty Hamiltonian - i.e. in this case exactly 
half of the stabilizer group [19] will anti-commute with 
each error and Wj = 2 N ~ n, ~ 1 Vj, where N is the num- 
ber of physical qubits and n; is the number of logical 
qubts. Alternatively one could scale the energy penalty 
to be exponential in the number of qubits: log a ~ O(N). 
However, both approaches require exponential resources 
because exponential energy must be added to the Hamil- 
tonian (either via the penalty a or via the number of 
stabilizer terms added). It is an interesting question as 
to whether such exponential suppression is possible with- 
out this concomitant cost. For this model of noise (weak 
and not strongly non-Markovian) we conjecture that this 
is not possible. 

Description of the dynamics as a classical rate 
equation is only possible when -Haqc = 0. When 
this is not the case, consider the traces in Eq. (4): 
Tr{E j PE j (t,T)Qp(t)Q} is still equal to Tr{QipQi}, 
but we can only bound the second trace as 
TT{P-E,-3j-(t,T)P/3(i)P} < P c (t). Although a rate 
equation for the populations is not possible in this case, 
it is clear that the mechanism by which the control 
suppresses leakage from the codespace remains the same. 
We would ideally like dP c (£)/dt = and since the trace 
quantities will not be zero in general, the alternative is 
to suppress the magnitude of the integral coefficients. 
This is exactly what the control Hamiltonian does by 
adding oscillatory components to the integrand defining 
the population transition rates, although the expressions 
for the suppressed leakage rates are not as simple as in 
the Haqc = case. 

Discussion - A practical consequence of the equivalence 
of the DD and EGP error suppression techniques is that 
it implies that DD can be used to emulate EGP. This 
is significant since EGP typically requires the addition 
of many-body terms (terms with a non-identity Pauli on 
more than two qubits of the system) to the system Hamil- 
tonian, which is usually impractical. In addition, as il- 
lustrated above the equivalence also prompts generaliza- 
tions of EGP where the penalty Hamiltonian coefficients 
(a) are not uniform and time dependent. Finally, our 
analysis of the limitations of error suppression in AQC 
motivates schemes for error correction within this model 
of quantum computation, while we have also shown that 



5 



implementations of error correction are non-trivial due to 
the always-on adiabatic evolution. In a follow-up paper, 
motivated by these observations, we will formulate meth- 
ods of error correction that are compatible with AQC. 
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SUPPLEMENTARY INFORMATION 



A. Appendix: Detectable errors commute with exactly half of the stabilizer group 

Theorem: All detectable error operators of an error detecting code must anti-commute with exactly half of the full 
stabilizer group and commute with the other half. 

Proof: The stabilizer group, Q, of an error detecting code is an abclian group generated by a set of stabilizer generators, 
S. The set of detectable errors, £ , is composed of all operators which anti-commute with at least one of the stabilizer 
generators. Because both the stabilizer group and the set of correctable errors are contained within the 7V-qubit Pauli 
group, all such operators either commute or anti-commute with one another. 

Any given detectable error operator, Ej, will partition Q into commuting and anti-commuting subsets Qj C G = 

{g G G\Ejg = ±gEj}. Because Ej is detectable, we can choose a generator, Sq, with which it anti-commutes. This 
generator defines an involutive map between the sets Gf 1 by g G Gf S° Q g G Gj ■ That is, if g commutes with Ej, 
then Sjg anticommutes with Ej, and visa- versa. The sets Gf are then isomorphic, so are the same size. 



B. Appendix: Derivation of encoded AQC master equation 

In this Appendix we outline the derivation of a master equation describing the effective adiabatic evolution when the 
qubits are coupled to uncontrolled (environmental) degrees of freedom. By employing fewer approximations than in 
the derivation of the conventional Lindblad master equation this reduced dynamics is able to capture the modification 
of system-environment coupling, and hence decoherence, by controls such as dynamical decoupling or energy penalty 
terms. 

We begin with a Hamiltonian acting on the combined Hilbert space of encoded system and environment (H sy s ( ^'Hcnv) 
of the form: 

H(t) = H AQC (t) + H c (t) +J2 E i® Bj + H B (B.l) 

i=i 

where all these quantities are defined in the main text (Eq. (1)). Define an interaction picture with respect to HAQc(t), 
H c {t) and H B as: A(t) = W(t,0) AU(t,0), where 

W(t,0) = e+ /o< dsR ^ c{s)+Hc{s)+HB (B.2) 

and the subscript + indicated positive time ordering of the expocntial. A particularily important property of these 
Hamiltonians, which we will utilize later, is that they all commute. That is, [Haqc(s), Hc(s')] = Vs, s', because 
Haqc only contains logical operators and He only contains stabilizer terms. And obviously Hb commutes with 
the other two terms. This property implies that this interaction picture transformation factors into: U(t, 0) = 

Y[ n U n {t, 0) = n„ e+ f ° dsHn(s) with n G {AQC, C, B}. 

Let q be the combined density matrix of system and environment, i.e. is a normalized trace-class operator in 
H sys ®^cnv A formal solution to the von-Neumann equation in the interaction picture, = —i[H(t), g{t)], is [18]: 



dg(t) 
dt 



= -i[Hj(t), ~g(0)} - / [Hi(t), [H!(s), Q(s)]}ds (B.3) 



where Hj = Ylj^i-^j ® Bj. We will assume that the weak system-environment coupling does not perturb the 
environment from its equilibrium state at timescales that we resolve, and hence g(s) ss p(s) (g> a eq , a tensor product 
of the system density matrix p(s) = Tr env {g(s)}, and the environmental equilibrium density matrix. This allows us 
to derive a time-convolution master equation for the system density matrix [18]: 

dt ~ LTcnvi dt 1 

Ne ^ f 

= -J2 ds[C k3 (t,s)E k (t)E J (s)p(s)-C jk (s,t)E k {t)p(s)E J {s) 

- C k j(t,s)Ej(s)p(s)E k (t) + Cj k (s,t)p(s)Ej(s)E k (t)^ (B.4) 
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with Ckj(t,s) = Tr onv {Bk(t)Bj(s)(7 e q} is the quantum correlation function of the environment. To obtain this 
expression we have assumed that Tr env {Bj(t)a eq } = Vj - i.e. the average interaction force on the bath equilibrium 
state is zero. We assume that the environment is stationary, implying that this correlation function is only dependent 
on the time difference t = t — s. This simplifies the master equation to: 



dm 
dt 



= - E / dT ( C K ( T ^k (t)E, (t - r)p(t - t) - C* kj (r)E k (t)p(t - r)E 3 (t - r) 
„• i.-i Jo \ 

- C kj (r)E 3 (t - r)p(t - r)E k (t) + C* kj (r)p(t - r)E 3 (t - r)E k (t) 



3,k=l ■ 



(B.5) 



The final approximation we make is typically referred to as the first Markov approximation and replaces p(t — r) with 
p(t) in the integrals above [18]. This amounts to assuming that the change in the system state (in the interaction 
picture, and therefore due to the weak system-environment coupling) is negligible on the timescale set by the decay 
of the environment correlation function. Therefore this formalism is valid for moderately fast-relaxing or weakly- 
coupled environments. Finally, we will restrict out analysis to uncorrelated environments for the system qubits, that 
is C k j(r) = 5 k jCj(r). The analysis that follows can be generalized to correlated environments but we will not do so 
here. 

We rewrite this resulting master equation in an interaction picture with respect to the control Hamiltonian only. In 
the dynamical decoupling literature this is known as the toggling frame, and is defined by: A(t) = W^»(t, 0) A Uc(t, 0). 
The transformation required to move into this frame is particularly easy in this case because as noted above the 
stabilizer properties result in the factoring of the full interaction picture transformation unitary. In the toggling 
frame: 



dp(t) 
dt 



-i[HAQ C (t),p{t)] 



£/ dr( Cj {t)E 3 (t)Ej (t, r)p(t) - C* {t)E 3 {t)p(t)% (t, r) - C 3 {r)% (t, r)p{t)E 3 (t) + C* (r)/J(t)* ■ (t, r)E 3 (t) 



^AQc(M 



where Sj(t, r) = U A Qc(t, t-r) E 3 (t 

Now consider the change in the codespace (the no-error subspace) population as a result of these dynamics. As in 
the main text we define P as the projector onto the codespace, and Q = I — P. Then the change in the codespace 
population is = Tr{P^|P}. To evaluate this quantity, we will first insert identities in the form P + Q around 
p(t), resulting in: 

dP c (t) 
dt 



-C*{t) 



PE 3 (t)E 3 (t, r)Pp(i)P + PE 3 (t)E 3 (t, r)Qp(i)P 



PE j (t)Qp(t)QE j (t,T)P 



PE 3 (t,T)Qp(t)Q E j( t )' p 



+C?(r) 



Pp(t)PE 3 (t, t)E 3 (t)P + Pp(t)QE 3 (t, t)E 3 (t)P 



(B.6) 



The first of these is by definition 
i.e. Haqc(s) and Hc{s) cannot 



where we have used the identities: PQ = 0, PEjP = Vj, and PE 3 (t, r)P = Vj. 
and the others follow from the properties of the Hamiltonian and error operators - 
move states between the subspaces defined by P and Q, and E 3 applied to any state in P results in a state in Q. 
The term Tr {PEj(t)Ej(t, r)Qp(t)P} and its conjugate also evaluate to zero although it is slightly more involved to 



see why. The reason is that E 3 (t)E 3 (t, r) = W c {t, Q)E 3 U c {t 1 0)£Yaqc(*) t 



r,0)E j U c (t 



T,0)U{ QC (t,t 



contains two applications of E 3 interleaved with unitary evolution that does not connect different stabilizer sectors 
and hence cannot connect P and Q subspaces. Hence, this master equation simplifies to: 



dPcjt) 
dt 



-C*{r) 



PE 3 (t)E 3 (t, T )Pp{t)V 



+ C*Ar) 



Pp{t)PE 3 {t,r)E 3 {t)P 



PE 3 (t)Qp(t)QZ 3 (t,T)P 



PE 3 (t,T)Qp(t)QE 3 (t)P 



(B.7) 
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At this point wc employ a critical property of the control Hamiltonian: that it modulates the system-environment 
interaction. Using the expressions for toggling frame error operators in Eq. (2) and Eq. (3) allows us to simplify the 
equation of motion for codcspacc population to: 



dP c (t) 
dt 



2 YlJ o <^Re jCj(r)m j (t,T)Tr 



PS j (t,T)Q/J(t)QB i P 



P^H,-(i,r)Pp(i)P 



(B. 



where S(i, r) = Waqc(^ — T )^j ^aqc(^ — t )- I n * n ^ s equation the effects of control (error suppression) are 
encapsulated in the modulation functions rrij(t, r). For the two error suppression techniques, as shown in the main 
text, these functions take the form: 



< GP (i,r) 



,DD 



(t,r) = (-lJfW-fC*-'-) 



(B.9) 



where u>j is the number of stabilizer terms in the EGP penalty Hamiltonian that anti-commute with the error Ej. 
p(t) is the DD coefficient defined in the main text. To write the modulation function for EGP we have exploited the 
property 



S m ,-B,} = " 



(B.10) 



which follows from the fact that all states in the codespace are eigenvalue +1 eigenstates of the stabilizers. 

Since the correlation function decays with t, the value of the integrand at small values of r are the most important. 
And if we assume that Haqc varies slowly with respect to time, we can approximate 

UA Q c(t,t- r) = e-'A-^AQoM ^ e - iT H AQC (t) (b.ii) 

and hence approximate H(i,r) ~ E(t,r) = e~ lrHAQC ^ Ej e lTH ^Qc(t)_ Using this results in the final form of the 
population master equation used in the main text. 

The modulation functions given in Eq. (B.9) display some degree of asymmetry between the EGP and DD error 
suppression techniques because while m DD depends on times t and r, m p only depends on time r. We stress that 
this is only because we have restricted ourselves to the case of constant, uniform energy penalty a. As mentioned in 
the main text, a more general formulation of EGP (from which a unification with DD is even more straightforward) 
would allow for a to be time dependent: HQ GP (t) = — J3 m =i a j^)^m- In this case, 

mJ GP (t,T) = (e«)xW-x(*-r) (B.12) 

with x(t) = £{s E }=o Jo ^ s a *™( s )- I n thi s more general formulation the similarity between DD and EGP is even 
more evident. 



C. Appendix: Full rate calculation for Ohmic spectral density 



The example environment considered in the main text is a classical bath with exponentially decaying correlation. 
Here we generalize this to a true quantum environment and explicitly demonstrate that the controlled suppression 
of population leakage from the codespace holds in this case too. Consider a damped harmonic environment with an 
Ohmic spectral density with Lorentz-Drude regularization: J(ui) = 2Er^uj/ (w 2 + 7 2 ), where Er is the "reorganization 
energy" which quantifies the total system-environment coupling strength, and 7 is the inverse of the environment 
correlation timescale. The quantum correlation function for an environment with such a spectral density is: 



1 \ _, t ^4£ R7 



where /3 = 1/fcgT is the inverse temperature and v K = =52 are the Matsubara frequencies. For fast decaying 
correlations (large 7) the terms in the summand decay quickly and it is customary to truncate the sum at finite 
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k. Assuming that the error suppression technique is EGP and computing the rates in the population master/rate 
equation yields: 

ao[bo — je^ 7t cos(2aWjt — ^-) + 2awje~~ /t sin(2awjt — 4^)] 
Tj W = {2aw ] ) 2 + 1 2 

E= a K [b K — i> K e~ VKt cos(2aujjt) + 2aWje~ VKt svn{2awjt)} . . 

and similarly for rj (t). Here a%, hi are irrelevant constants. As in the case of a classical model of the environment, the 
suppression of these transition rates is achieved by two mechanisms: (i) the suppression term 2awj in the denominator 
decreases the overall rate of population leakage, (ii) the same term increases the oscillation frequency of the sinusoidal 
functions in the numerator, thus decreasing the magnitude of integrals of rf{t). 



